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Nonequilibrium kinetic Ising models evolving under the competing effect of spin flips at zero 
temperature and Kawasaki-type spin-exchange kinetics at infinite temperature T are investigated 
here in one dimension from the point of view of phase transition and critical behaviour. Branching 
annihilating random walks with an even number of offspring (on the part of the ferromagnetic 
domain boundaries), is a decisive process in forming the steady state of the system for a range of 
parameters, in the family of models considered. A wide variety of quantities characterize the critical 
behaviour of the system. Results of computer simulations and of a generalized mean field theory are 
presented and discussed. 
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I. INTRODUCTION 



The Ising model is a well known static, equilibrium model. Its dynamical generalizations, the kinetic Ising models, 
were originally intended to study relaxational processes near equilibrium states |]l],D. Glauber introduced the single 
spin-flip kinetic Ising model, while Kawasaki constructed a spin-exchange version for studying the case of conserved 
magnetization. To check ideas of dynamic critical phenomena |^], kinetic Ising models were simple enough for pro- 
ducing analytical and numerical results, especially in one dimension, where also exact solutions could be obtained, 
and thus have proven to be a useful ground for testing ideas and theories [Q. On the other hand, when attention 
turned to non-equilibrium processes, steady states, and phase transitions, kinetic Ising models were again near at 
hand. Nonequilibrium kinetic Ising models, in which the steady state is produced by kinetic processes in connection 
with heat baths at different temperatures, have been widely investigated and results have shown that various 

phase transitions are possible under nonequilibrium conditions, even in one dimension (Id) (for a review see the article 
by Racz in Ref. §). 

Also in the field of ordering kinetics, universality is strongly influenced by the conserved or non-conserved char- 
acter of the order parameter. The scaling behaviour and characteristic exponents of domain growth are of central 
importance. Kinetic Ising models offer a useful laboratory again to explore the different factors which influence these 
properties JlOf . 

Most of these studies, however, have been concerned with the effects the nonequilibrium nature of the dynamics 
might exert on phase transitions driven by temperature. In some examples nonlocal dynamics can generate long-range 
effective interactions which lead to mean-field (MF)-type phase transitions (see, e.g., Ref. §). 

A different line of investigating nonequilibrium phase transitions has been via branching annihilating random walk 
(BARW) processes . Here particles chosen at random carry out random walks (with probability p) and annihilate 
pairwise on encounter. The increase of particles in ensured through production of n offspring, with probability 1 — p. 
Numerical studies have led to the conclusion that a general feature of BARW is a transition as a function of p 
between an active stationary state with non-zero particle density and an absorbing, inactive one in which all particles 
are extinct. The parity conservation of particles is decisive in determining the universality class of the phase transition. 
In the odd-n case the universality class is that of directed percolation (DP) while the critical behavior is different 
for n even . The first numerical investigation of this model was done by Takayasu and Tretyakov . This new 
universality class is often called parity conserving (PC) (however, other authors call it the DI [Q or BAWe class p^); 
we will also refer to it by this name. A coherent picture of this scenario is provided from a renormalization point of 
view in [l^ . 

The first example of the PC transition was reported by Grassberger et al. who studied probabilistic cellular 

automata. These Id models involve the processes k ^ 3k and 2k ^ {k stands for kink), very similar to BARW with 
n — 2. The two-component interacting monomer-dimer model introduced by Kim and Park p[ represents a more 
complex system with a PC type transition. Other representatives of this class are the three-species monomer-monomer 
models of [^, and a generalized Domany-Kinzel SCA which has two absorbing phases and an active one. 

A class of general nonequilibrium kinetic Ising models (NEKIM) with combined spin fiip dynamics at T = 
and Kawasaki spin exchange dynamics at T = oo has been proposed by one of the authors [ pO| in which, for a 
range of parameters of the model, a PC-type transition takes place. This model has turned out to be very rich in 
several respects. As compared to the Grassberger CA models, in NEKIM the rates of random walk, annihilation 
and kink-production can be controlled independently. It also offers the simplest possibility of studying the effect 
of a transition which occurs on the level of kinks upon the behaviour of the underlying spin system. The present 
review is intended to give a summary of the results obtained via computer simulations of the critical properties of 
NEKIM Dynamical scaling and a generalized mean field approximation scheme have been applied in the 

interpretation. Some new variants of NEKIM are also presented here, together with new high-precision data for some 
of the critical exponents. 



II. NON-EQUILIBRIUM KINETIC ISING MODEL (NEKIM) 
A. Id Ising model 

Before going into the details of NEKIM, let us recall some well known facts about the one-dimensional Ising 
model. It is defined on a chain of length L: the Hamiltonian in the presence of a magnetic field H is given by 
Ti. = 3D — JJ2i ■SiSj+i ~ H Si, {si = ±1). The model does not have a phase transition in the usual sense, but when 
the temperature T approaches zero a critical behaviour is shown. Specifically, pt = plays the role of in Id, 

and in the vicinity of T = critical exponents can be defined as powers of pr- Thus, e.g., the critical exponent v of 
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the coherence length is given through ^ oc pt~'^ ■ Concerning statics, from exact solutions, keeping the leading-order 
terms for T — > 0, the critical exponents of {ksT times) the susceptibility, the coherence length and the magnetisation 
are known to be 7^ = = 1/2, and (3^ — 0, respectively. We can say that the transition is of first order. Fisher's 
static scaling law js — dv — 2(3^ is valid. 

The Ising model possesses no intrinsic dynamics; with the simplest version of Glauber kinetics (see later) it is 
exactly solvable and the dynamic critical exponent Z defined via the relaxation time of the magnetization as r = tq^'^ 
has the value Z = 2. (Here tq is a characteristic time of order unity.) 

Turning back to statics, if the magnetic field differs from zero the magnetization is also known exactly; at T = 
the solution gives: 

m{T = 0,H) = sgn{H). (1) 
Moreover, for ^ ^ 1 and H/kT <^ 1 the exact solution reduces to 

TO ~ 2/1^ ; /i = H/keT. (2) 

In scaling form one writes: 

m^r'^gihi^) (3) 

where A is the static magnetic critical exponent. Comparison of eqs. (^ and (||) results in (3^ — and A — v . It 
is clear that the transition is discontinuous at if = as well [i.e., upon changing the sign of H in Eq. (0)]. In the 
following the order of limits will always be meant as: (1) H ^ 0, and then (2) T ^ 0. 



B. Definition of tlie model 



A general form of the Glauber spin-flip transition rate in one-dimension for spin Si sitting at site 
I is [| (s. = ±l): 

= ^(1 + (Ss,_is,+i) (1 - ^s,(s,_i + s,+i)) (4) 

where 7 tanh 2 J/fcT , with J denoting the coupling constant in the ferromagnetic Ising Hamilto- 
nian, r and S are further parameters which can, in general, also depend on temperature. (Usually 
the Glauber model is understood as the special case (5 = 0, r = 1, while the Metropolis model 
is obtained by choosing r = 3/2, 6 = -1/3). There are three independent rates: 

r 

Wsame = 77(1 + ^ 7) (5) 
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r 



Woppo = + + (6) 



Windtf = ^{l~S), (7) 

where the subscripts same, etc., indicate the three possible neighborhoods of a given spin(tTT,iTi 
and tti, respectively). In the following T = will be taken, thus 7 = 1, Wsame = and r, S are 
parameters to be varied. 
The Kawasaki spin-exchange transition rate of neighbouring spins is: 

Wi,i+i{si, Si+i) = ^(1 - SiSi+i)[l - ^(si_iSi + (8) 
At T = 00 (7 = 0) the above exchange is simply an unconditional nearest neighbor exchange: 

Wi^i+l = ^Pex[l - SiSi+l] (9) 
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where pex is the probabihty of spin exchange. 

The transition probabihties in Eqs. (^) and are responsible for the basic elementary processes 
of kinks. Kinks separating two ferromagnetically ordered domains can carry out random walks 
with probability 

OC 2w7i„dj/ = r(l - (5) (10) 

while two kinks at neighbouring positions will annihilate with probability 

Pan OC Woppo r(l + S) (11) 

{wsame Is respouslble for creation of kink pairs inside of ordered domains at T ^ o). 

In case of the spin exchanges, which also act only at domain boundaries, the process of main 
importance here is that a kink can produce two offspring at the next time step with probability 

Pk->3k OC Pex- (12) 

The above-mentioned three processes compete, and it depends on the values of the parameters r, 
6 and pex what the result of this competition will be. It is important to realize that the process 
k ^ 3k can develop into propagating offspring production only if pmi > Pan, i.e., the new kinks are 
able to travel on the average some lattice points away from their place of birth, and can thus avoid 
immediate annihilation. It is seen from the above definitions that (5 < is necessary for this to 
happen. In the opposite case the only effect of the k ^ 3k process on the usual Ising kinetics is 
to soften domain walls. This heuristic argument is supported by simulations as well as theoretical 
considerations, as will be discussed in the next section. 

C. Phase boundary of the PC transition in NEKIM 

In all of our investigations we have considered a simplified version of the above model, keeping 
only two parameters instead of three by imposing the normalization condition prw + Pan + Pk~-,3k = l, 
i.e., 

2r + pe. = 1. (13) 

In the plane of parameters pex and S = phase diagram, obtained by computer simula- 

tion, is as shown in Fig. |l|. There is a line of second-order phase transitions, the order parameter 
being the density of kinks. The phase boundary was obtained by measuring p{t), the density of 
kinks, starting from a random initial distribution, and locating the phase transition points via 
p{t) OC with a = .28 ± .01. The lattice was typically L = 2000 sites; we generally averaged over 
2000 independent runs. The dotted vertical line in Fig. |l] indicates the critical point that we 
investigated in greater detail; further critical characteristics were determined along this line. 
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FIG. 1. Phase diagram of the two-parameter modeL The active phase hes above the boundary. The dotted vertical hne 
indicates the critical point that we investigated in greater detail. 



The line of phase transitions separates two kinds of steady states reachable by the system for 
large times: in the Ising phase, supposing that an even number of kinks are present in the initial 
states, the system orders in one of the possible ferromagnetic states of all spins up or all spins 
down, while the active phase is disordered, from the point of view of the underlying spins, due to 
the steadily growing number of kinks. 

For determining the phase boundary, as well as other exponents discussed below, we used the 
scaling considerations of Grassberger The density p{x,t,e) of kinks was assumed to follow 



the scaling form 

p{x, t, e) oc t-''<i>{ex^'''^ , ei^/"" ) (14) 

Here e measures the deviation from the critical probability at which the branching transition 
occurs, and are exponents of coherence lengths in space and time directions, respectively. 
By definition, z^n = u±Z, where Z is the dynamic critical exponent. <j){a,b) is analytic near a = and 
6 = 0. Using Eq. ( |I4D the following relations can be deduced. When starting from a random initial 
state the exponent f3 characterizes the growth of the average kink density in the active phase: 

p{e) = lim p{x,t,e) a (15) 

t — ^oc 

while the decrease of density at the critical point is given by 

p{t)oQt-°'. (16) 

The exponents are connected by the scaling law: 

P = v\f^a — v^Za (17) 

The phase boundary shown in Fig ^ was identified using Eq. (|T6| ) ||20|| . The initial state was 
random, with zero average magnetization. At the critical point marked on the phase boundary we 
recently obtained the value a = .280(5). This point was chosen in a region of the parameters {pex,5), 
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where the effect of transients in time-dependent simulations is small. Critical exponents were 
measured around the point: r = 0.35, 5c = -0.395 ±.001 and pex = 0.3. The deviation from the critical 
point e was taken in the 5-direction. It is worth noting that in all of our simulations of the NEKIM 
rule, spin-flips and spin-exchanges alternated at each time step. Spin-flips were implemented using 
two-sublattice updating, while L exchange attempts were counted as one time-step for exchange 
updating. The numerical value of the phase transition point is sensitive to the manner of updating 
(thus, e.g., with sequential updating in both processes, | 5c \ decreases by more than ten percent 
at the same values of r and Pex)- 

Concerning /3, we carried out quite recently a high-precision numerical measurement; the result 
is shown in Fig. |^. To see the corrections to scaling as well, we determined the effective exponent 
/3e//(e), which is deflned as the local slope of p{e) in a log- log representation between the data points 

where ei = 5c- 5i, providing an estimate 

/3=lim%/(e). (19) 

As shown in Fig. ^ the effective exponent increases slowly and tends towards the expected PC 
value as e 0. A simple linear extrapolation yields the estimate (3 = 0.95(2). 



1.1 




FIG. 2. Effective exponent /3e// near the critical point. Simulations were performed on a Id NEKIM ring of size L — 24000. 
Averaging was done after the steady state was reached for one thousand data points and 10^ independent samples. 



The transients mentioned above in connection with the time-dependent simulations show up 
for small times and are particularly discernible at the two ends of the phase boundary, where 
the spin-flip and spin-exchange processes have very different time scales. For this reason, near 
Pex = 1.0, where we get close to (5 = (the Glauber case), it is very difficult to get reliable numerical 
estimates, even for a; the exponent has been found to grow slowly with time, from a value near 
zero at small times. It is important to notice that according to Eq. (^3]), pcx = 1 is approached 
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together with r 0; thus ^ oo, as otherwise 5c = cannot be reached. Several runs with 
different parameter values [including cases without the restriction, Eq. (^3[)], have been performed 
in this Glauber limit; all show that the 6 = case remains Ising-like for all values of r and pex'- the 
exponent a tends to 0.5 for late times. We therefore conjectured that the asymptote of the 
phase boundary is J = for = 1.0, and thus that > Pan is a necessary condition for a PC-type 
phase transition. 

In a recent paper, using exact methods, Mussawisade, Santos and Schiitz confirmed the 
above conjecture. These authors study a one-dimensional BARW model with an even number of 
nearest- neighbour offspring. They start with spin kinetics; in their notation D, A and a are the rates 
of spin diffusion, annihilation and spin-exchange, respectively, and transform to a corresponding 
particle (kink) kinetics to arrive at a BARW model. (In their notation D = pb.w, A = pan and 

a = 2pex-) 

Using standard field-theoretic techniques, Mussawisade et al. derived duality relations 

A = A, (20) 
b = \ + a, (21) 
a = D-\ (22) 

which map the phase diagram onto itself in a nontrivial way and divide the parameter space into 
two distinct regions separated by the self-dual line : 

D = \ + a. (23) 

The regions mapped onto each other have, of course, the same physical properties. iN particular, 
the line a = maps onto the line D = X and the fast- diffusion limit to the limit a ^ oo. In 
the fast diffusion limit the authors find that the system undergoes a mean-field transition, but 
with particle-number fiuctuations that deviate from those given by the MF approximation. So the 
nature of the phase transition in this limit is not PC but MF. The observation that for large D, any 
small a brings the system into the active phase, together with duality, predicts a phase transition 
at = A, in the limit a ^ oo. The exact result of confirms the conjecture - stemming from 



our simulations - on the location of the phase transition point in the phase diagram of NEKIM 
in the limit pex 1, which is indeed at 6 = 0. The dual limit D ^ oo covers the neighbourhood 
of the point S = p^x = 0. This result predicts an infinite slope of the phase boundary in this 
representation, at this point, as is also apparent in Fig. ^ It is important to notice that the duality 
transformation does not preserve the normalization condition of NEKIM, eq([ll), used in all of our 



numerical studies. The self-dual line of |27] D = A + a, which in terms of our parametrization reads 



6 = (24) 

1 -Pex 

does not conserve the above-mentioned normalization condition either. In the phase diagram of 
Fig. ID it corresponds to a line which starts at (0,0) and ends at (1/3,-1) in the {pex,S) plane. 

Finally, in ref. a proportionality relation between the (time-dependent) kink density for a 
half-filled random initial state, and the survival probability of two single particles in an initially 
otherwise empty system has also been derived. Concerning the respective critical exponents a and 
S (the latter being defined through the critical behaviour of the survival probability P{t) oc t^^), 
they have been shown to coincide within the error of numerical simulations in , as required by 
scaling. 



7 



D. Long-range initial conditions 



The symmetry between the two extreme inital cases in time dependent simulations (single seed 
versus homogeneous random state) concerning the kink density decay inspired us to investigate 
initial conditions with long-range two-point correlations in the NEKIM model. This procedure 
was shown to cause continuously changing density decay exponents in case of DP transitions . 
In this way it has been possible to interpolate continuously between the two extreme initial cases. 

In the present case we have started from initial states containing an even number of kinks 
to ensure kink number conservation mod 2 and possessing kink-kink correlations of the form 
< kiki+x > a;~(^~'^) with cr lying in the interval (0, 1). Here h denotes a kink at site i. These states 
have been generated by the same serial algorithm as described and numerically tested in ref. . 

The kink density in surviving samples was measured in systems of L = 12000 sites, for up to 
tmax = 80000 time steps; we observe good quality of scaling over a time interval of three decades 
(80,80000) (see Fig. § [g^. 

As one can see, the kink density increases with exponent a ^ 0.28 for o- = 0, where in principle 
only one pair of kinks is placed on the lattice owing to the very short range correlations. Without 
the restriction to surviving samples we would have expected a constant, the above exponent 
corresponds to dividing the kink density by the survival probability,P(i) oc t^^ and the above value 
is in agreement with our former simulation results for 6 ||24|| . 

In the other extreme case for a = I the kink density decays with with an exponent a ~ -0.28, 



again in agreement with our expectations in case of a random initial state ||2^. In between the 
exponent a changes linearly as a function of a and changes its sign at cr = 0.5. This means that the 
state generated with cr = 1/2 is very near to the situation which can be reached for t ^ oo, i.e., the 
steady-state limit. 



'c 

03 

1^ 



0.1 



0.01 



0.001 




100 1000 10000 

t 



FIG. 3. log(pfci„fc(t)) versus log(f) in NEKIM simulations for initial conditions with a = 0, 0.1, 0.2..., 1 (from bottom to top). 

We mentioned at the end of the last section that the duality transformation of ref. predicts 
the equality of the above two extremal exponents (though it also follows from scaling). One can 

make the conjecture that in greater generality the duality transformation may connect also cases with initial 
two-point correlation exponents: a {1 — a). 



E. Variants of the NEKIM model 
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1. Varying exchange-range model 



Now we generalize the original NEKIM model by allowing the range of the spin-exchange to 
vary. Namely, Eq. (|^) is replaced by 

Wi^i+k = -jPex[l - SjSi+fc], (25) 

where i is a randomly chosen site and is allowed to exchange with Si+k with probability pex- Site 
k is again randomly chosen in the interval \ < R , R being thus the range of exchange. The spin-flip 
part of the model will be unchanged. We have carried out numerical studies with this generalized 
model ||21j| in order to locate the lines of Ising-to-active PC-type phase transitions. The method 
of updating was as before. It is worth mentioning, that besides k 3fc, the process k ^ 5k can 
also occur for i? > 3, and the new kink pairs are not necessarily neighbors. The character of the 
phase transition line at i? > 1 is similar to that for i? = 1, except that the active phase extends, 
asymptotically, down to Sc = 0. 

1.0 n 



ACTIVE PHASE 

\ 




pcx r7(i+r)' 

FIG. 4. Phase diagram of NEKIM for 5c{R,Pex) and a) 7i = 1, and b) J? = 3, as a function of pex (with 2r = 1 — Pex), and 
for c) r = 0.35 {pex = 0.3) as a function of {R/l + Rf (full line). 



This is illustrated in Fig. H, where besides i? = 1, the case i? = 3 is also depicted: the critical 
value of -5c is shown as a function of pex [curves a) and b)]. Moreover, -5c as a function of R is also 
shown at constant r = .35 . The abscissa has been suitably chosen to squeeze the whole (infinite) 
range of R between and 1 and for getting phase lines of comparable size (hence the power 4 of 
i?/(l + R) in case of curve c)). The phase boundaries were obtained by measuring p{t), the density 
of kinks, starting from a random initial distribution and locating the phase transition points by 
p{t) oc with a = 0.27 ± .04. We typically used a system oi L = 2000 sites and averaged over 500 
independent runs. 

Besides the critical exponents we have also determined the change of -5c with \/R numerically 
at fixed r = .35, pe^^ = .3 . Over the decade of i? = 4 - 40 we have found 

-5c « 2.0(l/i?)', (26) 
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which is reminiscent of a crossover type behaviour of equihbrium and non-equihbrium phase tran- 
sitions [Q, here with crossover exponent 1/2. It should be noted here, that to get closer to the 
expected 6cmf = 0, longer chains (we used I values up to 20000) and longer runs (here up to i = S-IO'^) 
would have been necessary. The former to ensure I > R and the latter to overcome the long 
transients present in the first few decades of time steps. In what follows we will always refer to 
the MF limit in connection with 1 (i-e., Pex/^ —> oo), for the sake of concreteness, but keep in 
mind that oo can play the same role. 



2. Cellular automaton version of NEKIM 



Another generalization is the probabilistic cellular automaton version of NEKIM, which consists 
in keeping the rates given in eqs.(0) and prescribing synchronous updating. Vichniac ||3^ investi- 
gated the question how well cellular automata simulate the Ising model with the conclusion that 
the unwanted feedback effects present at finite temperatures are absent at T = and the growth 
of domains is enhanced; equilibrium is reached quickly. In the language of Glauber kinetics, this 
refers to the S >0 case. 

With synchronous updating, it turns out that the T = Glauber spin-flip part of NEKIM itself 
leads to processes of the type k 3fc, and for certain values of parameter-pairs (r,S) with 6 < 0, 
a PC-type transition takes place. For r = .5 the critical value of (5 is 5c = -.425(5) , while e.g., for 
r = .35, (5c = -.535(5). For | (5 1 >| Jc I the phase is the active one,while in the opposite case it is 
absorbing. All the characteristic exponents can be checked with a much higher speed, than for the 
original NEKIM. The phase boundary of the NEKIM-CA in the (r, -S) plane is similar to that in 
Fig. ^except for the highest value of r = 1 Jc = cannot be reached; the limiting value is Sc = -.065. 
For the limit S = Vichniac's considerations apply, i.e., it is an Ising phase. The limiting case 
r ^ with Sc -1 is very hard from the point of view of simulations, moreover the case r = is 
pathological: the initial spin-distribution freezes in. The corresponding phase diagram is shown 
in Fig. I]. 
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FIG. 5. Phase diagram of the CA version of NEKIM . The phase transition points have been determined in the same way 
as for Fig.l 
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3. NEKIM in a magnetic field 



It has been shown earher, in the framework of a model different from NEKIM exhibiting a 
PC-type phase transition [|3l|], that upon breaking the Z2 symmetry of the absorbing phases, 
DP universality is recovered. In case of the Glauber model a similar situation arises with the 
introduction of a magnetic field. In the presence of a magnetic field the Glauber transition rates 
given in Sec. 2 are modified so: 



'^'indif 
^oppo 



ppo (1 hSi)^ 



h = tanh(— — ) 

KrT 



(27) 
(28) 

(29) 



Fig. ^ shows the phase diagram of NEKIM in the {h,5) plane ||2^, starting at the reference PC 
point for ft, = used in this paper {5c = -.395). 
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FIG. 6. Phase diagram of NEKIM in the {h, 5) plane in the presence of an external magnetic field. Parameters of the 
transition probabilities: F = .35,Pei = -3. Naturally, the phase diagram can be drawn symmetrically for negative values of h, 
as well. 



We have applied only random initial state simulations to find the points of the line of phase 
transitions (critical exponents: a = .17(2), j3 = .26(2), corresponding to the DP universality class, as 
expected.) It is seen, that with increasing field strength the critical point is shifted to more and 
more negative values of 5. 

Finally, it is worth mentioning a further possible varint of NEKIM - again without a magnetic 
field - in which the Kawasaki rate Eq. is considered at some finite temperature, instead of 
T = oo, but keeping T = in the Glauber dynamics. As lowering the temperature of the spin- 
exchange process acts against kink production, the result is that the active phase shrinks. For 
further details see Ref. [02 . 
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III. GENERALIZED MEAN-FIELD THEORY AND COHERENT ANOMALY EXTRAPOLATION 



Using the CA version of NEKIM described in the previous section we performed n-site cluster 
mean-field calculations (GMF) up to n = 6, and a coherent anomaly extrapolation (CAM), resulting 
in numerical estimates for critical exponents in fairly good agreement with those of other methods. 
With this method we could also handle a symmetry-breaking field that favors the spin-fiips in one 
direction. For this SCA model we applied the generalized mean-field technique introduced by 
Gutowitz and Dickman for nonequilibrium statistical systems. This is based on the 



calculation of transitions of n-site cluster probabilities -P„({si}) as 

where is a function depending on the update rules and the last equation refers to our solution in 
the steady state limit. At the fc-point level of approximation, correlations are neglected for n> k, 
that is, P„(si, s„) is expressed by using the Bayesian extension process [^,^ . 



WjZl''' Pk-l{si+j, ...,Sk-l+j 



= ; • (31) 



The number of independent variables grows more slowly than 2" - 1, owing to internal symmetries 
(refiection, translation) of the update rule and the 'block probability consistency' relations: 

P„(S1, S„) = ^ F„+i(si, Sn, Sn+l) , 
s„+i 

P„(S1, Sn) = ^ Pn+l{so, Si, Sn) • 
So 

The series of GMF approximations can now provide a basis for an extrapolation technique. One 
such a method is the coherent anomaly method (CAM) introduced by Suzuki [^^. According to 
the CAM (based on scaling) the GMF solution for a physical variable A (example the kink density) 
at a given level (n) of approximation - in the vicinity of the critical point, Sc - can be expressed 
as the product of the mean- field scaling law multiplied by the anomaly factor a„: 

A, =a„ (<5/<5,"-l)^«- , (32) 

and the n-th approximation of critical exponents of the true singular behavior (0„) can be obtained 
via the scaling behavior of anomaly factors: 

an ^ A^"^'■^^ (33) 

where we have used the A„ = {Sc/S^ - S^/Sc) invariant variable instead of e, that was introduced to 
make the CAM results independent of using 6^ or 1/Sc coupling ( [0). Here 6^ is the critical point 
estimated by extrapolation from the sequence approximations S^. Since we can solve the GMF 
equations for n < 6, we have taken into account corrections to scaling, and determined the true 
exponents with the nonlinear fitting form : 

where b and c are coefficients to be varied. 

As discussed in a previous section, by duality symmetry the pex = axis is mapped into the S = 
axis and the neighborhood of the {pex ^ l,S = o) tricritical point which can be well described by the 
mean-field approximations to the [pex = 0,6 = -l) point. Therefore this point is also a tricritical 
point and the mean-field results are applicable in the neighborhood. 
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A. Z-i symmetric case 



First we show the results of the GMF+CAM calculations without external field. For n = 1, we 
have the traditional mean-field equation for the magnetization: 



dm/dt = —5Tni{m^ — 1) 

and for the kink density : 

dp/dt = 26Vp{l-ip + 2p^) 
The stable, analytic solution of the kink density exhibits a jump at 5 = 0. 



pi(oo) = 



1/2 : 5 < 
(5 > 







The n = 2 approximation gives for the density of kinks, p{oo), the following expression: 



P2(oo) 



4 I' 



rw ^ Pan Prw Pan \J ^qPtw 2 Prw ^ Pan 2 Prw^ Pan Pan^ 2 Pj-w Pan^ 



2 (^2Pr'w'^ Prw Pan Pan^^ 



(35) 



(36) 



(37) 



for (5 < 0. For 5 >Q p{oo) = 0, i.e., GMF still predicts a first-order transition for 5 = 0; the jump size 
in p(oo) at (5 = 0, however, decreases monotonically with decreasing r, according to eq. (p7|). 





FIG. 7. p(c>o) as a function of 5 as obtained from GMF for n = 3, .., 6 in case of a: F = .35, b: F = .05 



The GMF equations can only be solved numerically for n > 2. We determined the solutions of 
the n = 3,4, 5,6 approximations for the kink density at i). F = 0.35 (Fig.^ a) and of the n = 3,4, 5 
approximations at ii). F = 0.05 (Fig. (0 b)). 

As we can see, the transition curves become continuous, with negative values for (5" (5" denotes the 
value of 5 in the n-th approximation for which the corresponding p{oo) becomes zero). Moreover, 
I (5" I increases with growing n values. As increasing n corresponds to decreasing mixing, i.e., 
decreasing pex, the tendency shown by the above results is correct. 

Figs. P a) and || b) show a quantitative - though only tentative - comparison between the results 
of GMF and the simulated NEKIM phase diagrams. The obtained GMF data for 5", corresponding 
to n = 3, 4, 6 (f = 0.35) are depicted in Fig. ^ a) as a function of l/{n - 3), together with results of 
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simulations. The correspondence between n and pex was chosen as the simplest conceivable one. 
(Note that (5c 7^ is obtained first for n = 4.) The simulated phase diagram was obtained without 
requiring the normalization condition, Eq. (|13|), at fixed r = 0.35. In this case, the = limit, of 
course, is not reached and a purely second-order phase transition line can be compared with the 
GMF results (for n values where the latter also predicts a second-order transition). 




0.0 0.0 

0,0 0,2 0.4 0.6 0,8 1.0 0.0 0,2 0.4 0.6 0.8 1.0 

pex '/(n-3) pex 1/(n-3) 



FIG. 8. Comparison between simulation results and GMF predictions. GMF results for 5c are plotted as a function of 
l/(n — 3), while simulation results for 5^ are depicted as a function of Pea: at constant F, with 7? = 3, F = 0.35 in case (a), and 
(b) with i? = 1 for F = 0.05. 



Simulations for i? = 3 were found to lead to | (5c | values low enough to fit GMF data. The 
(polynomial) extrapolation of GMF data to n oo (corresponding to pex = , i-e., plain spin- flip), 
also shown in Fig. ^, could be expected to approach (5 -1. That this is not case may be due to the 
circumstances that upon increasing n i) GMF starts here from a first-order MF phase transition, 
which ii) becomes second-order, and iii) GMF should end up at the dual point discussed in Sec. 
II. C, while this symmetry is not included in the applied approximation procedure. Fig. ^ (b) 
shows the n = 3, 4, 5 results for F = 0.05, which are compared now with the R = 1 simulation data. 

The results of our GMF approximation useful for a CAM analysis are the n = 4, 5, 6 data, while 
the n = 3 result can be taken to represent the lowest-order MF approximation (with (5c*^^ = O) for 
a continuous transition (no jump in ps). For Sc we use a polynomial extrapolation. 

In the n 3 approximation the exponent p = 1.0064, thus (3mf ~ 1- Consequently, as the table 
below shows, the anomaly factor does not depend on n. This means, according to Eq. (|33|), that 
the exponent is estimated to be equal to the MF value /3 ~ /3mf = 1- 



TABLE L CAM calculation resuhs 



n 






4 


2.49043 


0.01083 


5 


1.81022 


0.01074 


6 


1.45766 


0.01079 
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B. Broken Zi symmetry case 



The transition probabilities of NEKIM, in the presence of en external magnetic field H have been 
given in eqs.(27)-(29), in the previous section. Now we extend the method for the determination 
of the exponent of the order-parameter fluctuation as well : 

X(e) = L(< > - < p >^) ^ 6--^" . (38) 

The traditional mean-field solution [n = l) results in stable solutions for the magnetization : 

i/ (5 < and /iV<5^ < 1 (39) 

(40) 

(41) 
(42) 



mi 



mi = sgn(/i), otherwise. 

and for the kink-concentration : 

pi(cx)) = i(l - (1)2), z/ ,5<0 and h^l8^<\ 
pi(oo) — 0, otherwise. 

For n > 1 the solutions can again only be found numerically. Increasing the order of approxima- 
tion, the critical point estimates (5" shift to more negative values similarly to the H ^0 case. The 
hni„^oo <5"(/i) values have been determined with quadratic extrapolation for h = 0.01,0.05,0.08, and 
0.1. The resulting curves for p„((5) and x„((5) are shown in Figs. ^ (a) and (b), respectively, for the 
case of /i = 0.1 in different orders n of the GMF approximation. 



0.3 



0.2 



0.1 



-0.8 



-0.6 



-0.4 



-0.2 



0.2 



0.1 



-0.8 



o 

III I II 



-0.6 -0.4 



-0.2 



FIG. 9. The kink density (a) and (b) the second moment of the kink density in the neighbourhood of the critical point Sc{h) 
for h = 0.1. The curves from right to left correspond to n = 1, ...,6 (level of GMF calculation). The points were determined 
with a resolution of 10~^ in S, in order to be able to extract CAM anomaly coefficients. 



Naturally, these curves exhibit a mean-field type singularity at the critical point : 

Pn ^p„{S/Sl^~ir^'- (43) 

Xn~X„(V'5?-l)"^"", (44) 

with Pmf ^ and -jmf = -1- The results of the CAM extrapolation for various /i- values are shown 
in the table below : 
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TABLE II. CAM calculation results 



h 


0.0 


0.01 


0.05 


0.08 


0.1 


DP 


PC 


P 


1.0 


0.281 


0.270 


0.258 


0.285 


0.2767(4) 


0.94(1) 


7 




0.674 


0.428 


0.622 


0.551 


0.5438(13) 


0.00 
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As to the case /i = 0, we could not determine the exponent 7„ because the low level GMF calcu- 
lations resulted in discontinuous phase transition solutions - which we cannot use in the CAM 
extrapolation - and so we had too few data points to achieve a stable non-linear fitting. Higher- 
order GMF solutions would help, but that requires the solution of a set of nonlinear equations 
with more than 72 independent variables. This problem does not occur for h^O; the above results 
- being based on the full set of approximations (n = 1, ...,6) - are fairly stable. 



IV. DOMAIN GROWTH VERSUS CLUSTER GROWTH PROPERTIES; HYPERSCALING 

In the field of domain-growth kinetics it has long been accepted that the scaling exponent of 
L{t), the characteristic domain size, is equal to 1/2 if the order parameter is not conserved. In the 
case of a Id Ising spin chain of length L, the structure factor at the ferromagnetic Bragg peak is 



defined as: S{0, t) = L[< > - < m >\ 



L 



[3, 



±l). If the conditions of validity of scaling 



are fulfilled 38 then 



(45) 



where now d = I and x = 1/2. Another quantity usually considered [|3g] is the excess energy 
A£;(t) = E{t) - Et { Et is the internal energy of a monodomain sample at the temperature of 
quench), which in our case is proportional to the kink density p{t) = x < Si 5(1 ^ SiSt+i) >. 



p{t) 



Lit) 



oc t 



(46) 



with y = 1/2 in the Glauber- Ising case, expressing the well known dependence on time of annihi- 
lating random walks (y = a in the notation of section 2.). We also measured S{t) at the PC point 
and found power-law behaviour with an exponent x = .575(5) ^2 . 



The second characteristic growth length at the PC point is the cluster size, defined through the 
square-root of < i?^(t) xxt^. The latter is obtained by starting either from two neighbouring kink 
initial states (see e.g. [0), or from a single kink [ p!8| , |20[] . Both length-exponents are, however 
connected with Z, the dynamic critical exponent, since at the PC transition the only dominant 
length is the (time-dependent) correlation length. It was shown in |jT^ and ||3^] that | = ;| follows 
from scaling for one-kink and two-kink initial states, respectively. In our simulations we also found 
2^ that X = l/Z, within numerical error, at the PC point, and thus 



X = z/2 



(47) 



follows. Exponents of cluster growth are connected by a hyperscaling relation first established 



by Grassberger and de la Torre |^ for the directed percolation transition. It does not apply to 
the PC transition in the same form [12|, where the dependence on the initial state (one or two 
kinks) manifests itself in two cluster-growth quantities: the kink-number N{t,e) oc f'f{et^^) and the 
survival probability P(t, e) oc t~^g{et~) (This S has, of course, nothing to do with the parameter 
S of NEKIM). In systems with infinitely many absorbing states, where 6 and 77 depend on the 
density of particles in the initial configuration, a generalised scaling relation was found to replace 
the original relation, 2(5 + 77 = z/2, namely [p |: 



2(/3 +/3) 



I Z 



277 



(48) 
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which is also apphcable to the present situation. In eq(^) /?' is defined through linit^oc P{t, e) oc e'' 
and the scahng relation = S'iy±Z holds. A prime on an exponent indicates that it may depend on 
the initial state. For the PC transition the exponent z has proven (within error) to be independent 
of the initial configuration, z' ^ z. 

In case of a single kink initial state all samples survive thus S' = 0, and via the above scaling 
relation also = 0. Eq(|48|) then becomes: + 2ri' = z. For t oo, however, limt^oo N{t, e) cx has 
to hold ( the steady state reached can not depend on the initial state provided samples survive). 
Thus using the above scaling form for N{t,e), 77 = -^^^^ follows, and the hyperscaling law can be cast 
into the form: 



4/3 
ly I Z 



(49) 



Starting with a two-kink initial state, however, p' = (3 and ry' = holds [0. With these values 
Eq. ( PSP leads again to Eq. (^). As z = 2/Z, Eq. (^) involves only such quantities, which are 
defined also when starting with a random initial state. Using Eq. (|l3), Eq. (^9|) becomes 



2y^l/Z = x (50) 

In this way the factor of two between x and y appearing at the PC point between the exponents of 
a spin-bound and a kink-bound quantity (and which has been found also for a variety of critical 
exponents, see p2|) could be explained as following from (hyper) scaling. 



V. EFFECT OF THE PC TRANSITION ON THE PROPERTIES OF THE UNDERLYING SPIN SYSTEM 



Time-dependent simulations, finite-size scaling and the dynamic early-time MC method [43 



have been applied to investigate the behaviour of the Id spin system under the effect of the PC 
transition. We found |j2^ that within error of the numerical studies, the dynamic critical exponent 
of the kinks Z and that of the spins Zc agree, and have the value: Z ^ Zc = 1.75. Thus in comparison 
with the Glauber-Ising value Z = 2.0, Z decreases, as do 7 and ly. ^ = v ^ .444 (instead of 1/2). 
The PC point is the endpoint of a line of first-order phase transitions (by keeping pex and r fixed 
and changing 5 through negative values to (5c), where (3 = still holds, as does Fisher's scaling law: 
^ = dv - 2(3. To obtain these values the PC point has been approached from the direction of finite 
temperatures of the spin-fiip process, by varying px = e~'^-^/^'^ . 

The dynamical persistency exponent 9 and the exponent A characterising the two-time 
autocorrelation function of the total magnetization under nonequilibrium conditions have also 
been investigated numerically at the PC point. It has been found that the PC transition has a 
strong effect: the process becomes non-Mar kovian and the above exponents exhibit drastic changes 



as compared to the Glauber-Ising case [23]. These results together with critical exponents obtained 



earlier in [E2| are summarized in Table III. 







7 


V 


Z 


e 


A 


Glauber-Ising 





1/2 


1/2 


2 


1/4 


1 


PC 


.00(1) 


.444(2) 


.444(2) 


1.75(1) 


.67(1) 


1.50(2) 



TABLE in. Simulation data for static and dynamic critical exponents for NEKIM 
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We have also investigated the effect of critical fluctuations at the PC point on the spreading 
of spins, and found the analogue of compact directed percolation (CDP) to exist. Compact 
directed percolation is known to appear at the endpoint of the directed percolation critical line 

ii 



of the Domany-Kinzel cellular automaton in 1 + 1 dimensions 



Equivalently, such a transition 



occurs at zero temperature in a magnetic field h, upon changing the sign of /i, in the one-dimensional 
Glauber-Ising model, with well known exponents characterising spin-cluster growth |52]. 




FIG. 10. Phase diagram of NEKIM in the {h, —S) plane. The chosen PC-point is at 5 = —0.395. For smaller values of —5, 
in the Glauber-Ising regime, the vertical line connecting the PC and the nonequilibrium Glauber Ising (NGI) points (at h = 0) 
consists of all CDP points, with its characteristic critical exponents. The simulations around the PC transition were done for 
the interval < h < 0.1. The other NEKIM parameters, in the whole plane, are: F = 0.35 and pex = 0.3. 



For the spreading process of a single spin s(i) = 1 in the sea of downward spins(s(j) = -1 for 
j 7^ i), the exponents 6s, r/s and Zs are defined at the transition point by the power-laws governing 
the density of I's Us oc the survival probability Ps{t) oc t~^% and the mean-square distance of 
spreading < R'^{t) >^ V"'. We obtained = 0, (5 = 1/2 and 2 = 1, in accord with the hyperscaling 
relation in the form appropriate for compact clusters: 



T^s+5s^ Zs/^ (51) 

At the PC point we found that the characteristic exponents differ from those of the CDP tran- 
sition, as could be expected. But basic similarities remain. Thus the transition which takes place 
upon changing the sign of the magnetic field is of first-order and its exponents satisfy Eq. ([51|). 
Accordingly it can be termed as 'compact', and we call it the compact parity-conserving transition 
(CPC). For the phase diagram, see Fig.|T^. 

On the basis of Eq. , the ^-dependence of the magnetization in scaling form may be written 

as 

m{t,h) r^r^giht^) (52) 
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We have investigated the evolution of the nonequihbrium system from an almost perfectly mag- 
netized initial state (or rather an ensemble of such states) . This state is prepared in such a way 
that a single up-spin is placed in the sea of down-spins at Ljl. Using the language of kinks this 
corresponds to the usual initial state of two nearest neighbour kinks placed at the origin. At the 
critical point the density of spins is given as 

n,{tM)^f^''gx{hf^) (53) 

for the deviation of the spin density from its initial value, Us = m{t,h) - m(0). The results are 
summarized in Table IV. 







0s' 


A 


Vs 


Ss 


Zs 


NGI-CDP 





.99(2) 


1/2 


.0006(4) 


.500(5) 


1(= 2/Z) 


Cl'C 


.00(2) 


.i-Xi) 


.19(1) 


.288(1) 


.287(:!) 


1.11(=2/Z) 



TABLE IV. Spin-cluster critical exponents for NEKIM in a magnetic field The NGI-CDP data are results of simulations at 
the point 5 = 0, pex = 0.3 
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Concerning static exponents, only the magnetic exponent A remains unchanged under the effect 
of the critical fluctuations. This is not an obvious result (see, e.g., the coherence- length exponent 
above), (is , characterizing the level-off values of the survival probability of the spin clusters is 
deflned through 

YiTCL Ps{t,h) o^hf^^' ] (54) 
this is a new static exponent connected with others via the scaUng law 

= ^ (55) 
which is satisfled by the exponents obtained numerically, within error. 



VI. DAMAGE-SPREADING INVESTIGATIONS 



While damage spreading (DS) was introduced in biology [g^ it has become an interesting topic 



in physics as well [|i5|-pD|. The main question is if damage introduced in a dynamical system 
survives or disappears. To investigate this the usual technique is to construct a pair of initial 
conflgurations that differ at a single site, and let them evolve with the same dynamics and external 
noise. This method has been found very useful for accurate measurements of dynamical exponents 



of equilibrium systems [51 



In |^5| we investigated the DS properties of several Id models exhibiting a PC-class phase 
transition. The order parameter that we measured in simulations is the Hamming distance between 
the replicas: 



Dii) = (\il\^^-^\\) (56) 



=1 



where Si and the replica s\ may denote now spin or kink variables. 

If there is a phase transition point, the Hamming distance behaves in a power law manner at 
that point. In case of initial replicas differing at a single site (seed simulations) this looks like: 

D{t) oc f , (57) 
Similarly the survival probability of damage variables behaves as: 

Ps{t)(xt-^ (58) 
and the average mean square distance of damage spreading from the center scales as: 

R^it) oc . (59) 

In case of the NEKIM model we have investigated the DS on both the spin and kink levels. We 
found that the damage-spreading point coincides with the critical point, and that the kink-DS 
transition belongs to the PC class. The spin damage exhibits a discontinuous phase transition, 
with compact clusters and PC-like spreading exponents. The static exponents determined by 
flnite-size scaling are consistent with those of spins of the NEKIM model at the PC transition 
point and the generalised hyperscaling law is satisfled. 

By inspecting our large data set for various models exhibiting multiple absorbing states, we are 
led to a flnal conjecture: BAWe dynamics and Z2 symmetry of the absorbing states together form the necessary 
condition to have a PC- class transition. 
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VII. SUMMARY 



The present paper has been devoted to reviewing most of the authors' investigations of a nonequi- 
hbrium kinetic Ising model (NEKIM). 

NEKIM has proven to be a good testing ground for ideas, especially connected with dynamical 
scaling, in the field of nonequilibrium phase transitions. The field covered here belongs to the class 
of absorbing-state phase transitions which are usually treated on the level of particles, the best 
known examples being branching annihilating random walk models. The phase transition on the 
level of kinks in the Ising system belongs the the so-called parity-conserving universality class. The 
NEKIM model offers the possibility of revealing and clarifying features and properties which take 
place on the level of the underlying spin system as well. The introduction of a magnetic field in 
the Ising problem has also led to understand different features, e.g., the change of the universality 
class from PC to DP in the critical behavior at (and in the vicinity of) the phase transition of 
kinks. 

Investigations of various critical properties have been carried out mainly with the help of nu- 
merical simulations (time-dependent simulations, finite-size scaling and dynamic early-time MC 
methods have been applied). In addition, a generalized mean-field approximation was used, par- 
ticularly in the neighborhood of a limiting situation of the phase diagram of NEKIM, where the 
usually second-order phase transition line ends at a MF-type first-order point. Some of the critical 
exponents, especially that of the order parameter, could be predicted to high accuracy in this way. 

While from the side of numerics a lot of information has accumulated over the last five years 
concerning the PC universality class, further investigations are needed at the microscopic level. 
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